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Abstract. We will expose in this paper our advances towards a proof of 
the equivalence between FRW background expansion, during some period 
of time that contains primordial inflation, and the statistical isotropy of 
the primordial curvature perturbation ( at the end of this period of time. 
Our motivation rests on the growing interest in the existence of a preferred 
direction in the Universe hinted by the continuous presence of anomalies 
in the CMB data. 


1. Introduction 

Cosmology is based on a, once believed, sacred principle: the homogeneity and 
isotropy at large scales. The actual meaning of this statement is that perturba¬ 
tions in the energy density distribution at large scales, and in the CMB temper¬ 
ature, are statistically homogeneous and isotropic, i.e., their n-point correlators 
in real space are invariant under spatial translations and spatial rotations. How¬ 
ever, the CMB data releases have consistently presented a possible indication of 
the existence of a preferred direction in the Universe in relation with the different 
anomalies in the data (Ade et. ah, 2014; Bennett et. ah, 2013). Names as “the 
axis of evil” employed to describe some of the anomalies are evidence of both 
the idea of a preferred direction and the position of the cosmologists community 
about the violation of the sacred principle. There are two ways of thinking of a 
preferred direction in the Universe, one at the background level and the other at 
the perturbative level: at the background level we can have non-FRW expansion 
which can be associated, for instance, to the presence of a shear in the metric 
or an anisotropic curvature; whereas at the perturbative level we can have sta¬ 
tistical anisotropy in the perturbations, especially in the primordial curvature 
perturbation (). It is quite reasonable to think that a non-FRW background 
expansion always feeds the statistical anisotropy; moreover, sometimes a FRW 
background expansion is assumed in the analysis of the statistical anisotropy. 
Notwithstanding, there is no any proof yet regarding the actual relation between 
these two characteristics. It is our mission in this paper to expose our advances 
towards such a kind of proof. We will sketch the main argument of a proof 
of the equivalence between FRW background expansion, during some period of 
time that contains primordial inflation, and the statistical isotropy of the pri¬ 
mordial curvature perturbation ( at the end of this period of time, at the level of 
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the power spectrum, leaving the most delicate points for a deeper analysis and 
discussion in a future publication. 


2. Statistical homogeneity and isotropy 


Quantum mechanics only allows us to predict probabilities of different outcomes 
after an experiment in an ensemble of systems, in contrast to classical mechanics 
which does allow us to predict the exact outcome after an experiment in just 
one element of the ensemble. Since the underlying physical mechanism in the 
generation of cosmological perturbations is of quantum nature, the cosmologists 
are more interested in studying the statistical properties of a perturbation map, 
say the CMB map or the galaxy distribution map. The way of doing this is 
via the n-point correlators of the perturbations in real space. Let’s define a 
scalar cosmological perturbation /S(x) in real space and make a Fourier integral 
expansion 

/3(x) = j , (1) 

where /3(k) is the Fourier mode function of /3(x). The n-point correlators of /3(x) 
are averages over the ensemble of the products /3(xi)/3(x2).../3(x„) where xi, X 2 , 
..., x„ represent different points in spac4I|: 


(/?(xi)/3(x2).../3(Xn)) 


d^k 2 d?kn i(ki.xi+k2.x2+...+k„.x„) ^ 

(27r)^ (27r)3 (27r)^ 

x(/3(ki)/3(k2).../3(k„)). (2) 


Thus, the correlation functions in real space may be studied via the correlation 
functions in momentum space. Let’s see now the meaning of statistical homo¬ 
geneity and statistical isotropy. 

Statistical homogeneity: Of course the perturbation map is not homoge¬ 
neous (i.e., it is not invariant under spatial translations), but it may be that 
the probability distribution function governing /3(x) is, which is called statistical 
homogeneity. This means that the n-point correlators in real space are invariant 
under translations in space, i.e. 


(/3(xi -k d)/3(x2 d).../3(xn d)) = (/3(xi)/3(x2).../3(x„)), (3) 


where d is some vector in real space establishing the amount of spatial trans¬ 
lation. The only way of achieving this, in view of Eq. m, is expressing the 
argument in the exponential function inside the integral as the addition of sev¬ 
eral terms of the form /(xj — Xj), which in turn is possible (but it is not the only 
possibility) if the n-point correlators in momentum space are proportional to a 
Dirac delta function: 


(/3(ki)/3(k2).../3(k„)) = (27r)3h3(ki2,„„)M^(ki, k2,..., k„) . (4) 


^The ensemble average inside the integral is over the Fourier mode functions only since they are 
the stochastic variables. 
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In the previous expression, ki 2 ...n means ki + k 2 + ... + k^, and the function 
M,g(ki,k 2 , ...,k„) is called the (n — l)-spectrum. 

Statistical isotropy: Once statistical homogeneity has been secured, in the 
form of Eq. Q, we ask about the invariance under spatial rotations (i.e. isotropy). 
Of course again, the perturbation map is not isotropic, but it may be that the 
probability distribution function governing /3(x) is, which is called statistical 
isotropy. This means that the n-point correlators in real space are invariant 
under rotations in space, i.e. 

(/3(xi)/3(x2).../3(x„)) = (/3(xi)/3(x2).../3(x„)) , (5) 

where Xj = IZ Xi, IZ being a rotation operator. To satisfy the above requirement, 
the (n — l)-spectrum must satisfy the condition 

...jkn) = M/3(ki,k2, ...,k„), (6) 

where the tildes over the momenta represent as well a spatial rotation, parame¬ 
terized by TZ, in momentum space. This condition has more explicit consequences 
in the spectrum (1-spectrum) and the bispectrum (2-spectrum): 

M^(ki,k2) = P^(ki,k2) = P;3(A:), (7) 

A7y3(ki, k2, ks) = P^(ki, k2, ks) = i?,g(A:i,/c2, fea), (8) 

where in the first line /c = |ki| = |k 2 |, and in the second line ki = |kj|. Starting 
from the trispectrum (3-spectrum), the condition in Eq. Q about statistical 
isotropy in all the higher-order (n — l)-spectra cannot be reduced to similar 
conditions to the ones in Eqs. ([7|) and ([8]), so that the minimal way of param¬ 
eterizing the (n — l)-spectra (with n > 4) will always be in terms of all the n 
wavevectors. The scalar nature of /3(x) is very important since, if it were a vector 
or a tensor, there would not be a way to make the n-point correlators in real 
space invariant under spatial rotations. In those cases, we relax the meaning of 
statistical isotropy and establish that it is present if the (n — l)-spectra of the 
scalar pertubations that multiply the respective polarization vectors or tensors 
satisfy Eq. ([6]). 


3. The separate universe assumption and the 5N formalism 

The separate universe assumption: This assumption refers to the behaviour 
of the Universe after smoothing on a specified comoving scale k~^, during the 
super horizon era k <C aH. It states that the spatial gradients, at most of order 
k/a, are negligible, which actually means that the Universe at each comoving lo¬ 
cation behaves as if it were homogeneous (Lyth & Liddle, 2009). Each smoothed 
region about each comoving location is then regarded as a separate homogeneous 
universe. By virtue of this, the form of the equations for the dynamical quanti¬ 
ties in each separate universe is the same as for the unperturbed quantities. The 
separate universe assumption is a powerful tool for dealing with perturbations in 
the very early Universe and has become one of the most employed methodologies 
as alternative to the standard cosmological perturbation theory. 
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The 6N formalism: The 5N formalism (Dimopoulos et. ah, 2009; Lyth &: 
Liddle, 2009) provides a powerful method to evaluate the primordial curvature 
perturbation C(x, t) in terms of the perturbations of the fields a few Hubble times 
after horizon crossing (corresponding to a flat slicing), and the derivatives of 
the unperturbed number of e-foldings with respect to 

the unperturbed fields evaluated at t*. According to this formalism, once the 
separate universe approach has been invoked, and a comoving threading has 
been established, the value of in a uniform energy density hypersurface at the 
final time t is given by the perturbation in the time integral of the local volume 
expansion rate 9 along a curve starting at an initial flat hypersurface at the time 
U: 

C(x, t) = (5iV(x, t, ti) - (hA^(x, t, U)). (9) 

Here, the bracket notation means a ensemble average (which corresponds to a 
spatial average if there is statistical homogeneity). In many inflationary scenar¬ 
ios, the number N of e-foldings depends only on the values of the fields at t* so 
we can write the curvature perturbation as an expansion in the perturbations of 
the fields at this time. Assuming the presence of just one scalar field and one 
vector field (the generalization to more fields is straightforward), we have: 


Cix, t) = SN{4>{x.),Ai{x.),t) = N^5(p +Ni5Ai + ^N^^{S(j)f + N^i5(j)6Ai + 

Y—NijSAi6Aj , ( 10 ) 

where 

.r -dN d^N d‘^N 

= Q^2 ’ = QA.Q^ ’ ’ ( 11 ) 

are the derivatives with respect to the scalar (p and the spatial components of 
the vector field A. 


4. FRW background expansion implies statistical isotropy 

We will make the following assumptions: 

1. The action is such that the FRW metric is an attractor in the background. 

2. The fields involved are just scalar and/or vector fields. 

3. The background expansion is FRW type during the whole time spanning 
from the beginning of inflation to the time when the curvature perturbation 
C is evaluated. 

By invoking the separate universe assumption, the form of the equations 
for the dynamical quantities at each comoving location is the same as for the 
unperturbed quantities. Thus, the field perturbation equations in momentum 
space do not depend explicitly on k. They depend on time, thus on /c*, but do 
not depend on the direction of k. Their solutions are, therefore, independent 
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of the direction of k except for the set of initial conditions However, the 

field perturbations are evaluated in the flat slicing, so, taking into account the 
assumption [3l the whole perturbed metric in this slicing is actually FRW which 
is conformally equivalent to Minkowski. This has as a consequence that the set of 
initial conditions, when quantizing, can be written as {dJJ = 0 :^ 0 ^} where is 
the respective annihilation operator. Following the usual procedure to calculate 
the power spectrum of the field perturbations, this implies that the latter are 
statistically isotropic (and statistically homogeneous) at the level of the power 
spectrum (since the do not depend on the direction of k). This is valid for all 
the relevant cosmological scales since the time of horizon exit for each of them 
is contained in the period of time defined in the assumption |3l 

Now, if the background expansion is FRW type, the fluid that fills the 
universe is perfect; therefore, by invoking assumption [21 the held conhguration 
must be given by any number of scalar helds and/or multiple randomly oriented 
copies (just one is enough) of a triad of mutually orthogonal vector helds with 
the same norm. This must be accomplished during the whole period of time 
dehned in assumption mB Thus, by employing the 6N formalism, we obtain 
(Gomez & Rodriguez, 2013) 

Pc(ki) = priki)[l + 4(ki • Ni)2 + ■ N 2 )' + 5^(kl • N3)2] , (12) 

where 

Pf°(A:i) = {N^fPs4,ih) + (NlfPlih) + {NffPlih) + iNf)^Pl{h ), (13) 
and 

„ - Pi(h)] 

=- ffW) -’ ‘ ’ 

P+{k) being the parity even spectrum of the n-th vector held, PJ^ngi^) being the 
longitudinal spectrum of the re-th vector held, and 


N 


n 


N, 


N. 


(15) 


where N„ is the vector formed by the derivatives of N with respect to the each 
spatial component of the n-th vector held. However, because of the symmetries 
of the held conhguration, we have 


Ni 

= j 

(16) 

N 2 

= ±j, 

(17) 

N 3 

= 

(18) 

(Nl? 


(19) 

Pi 

II 

II 

(20) 

pi 

long 

_ p2 _ p3 

^long -^long ’ 

(21) 


^This has as a consequence that the action must be consistent with the isotropy during the 
referred period of time, e.g., the masses for the members of each triad must be equal (although 
not necessarily the same masses among triads). 
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and, therefore, 

9\= 9 }= 9l= 9(,- (22) 

Thus, 

Pc(ki) = Pf°(fci)[l + 5c[(ki-i)' + (ki-j)' + (ki-z)2]], (23) 

= Pf°(fei)[l + 5d, (24) 

where the last step is valid because of the director cosines property of a unit 
vector. We conclude then that the power spectrum of Q is actually indepen¬ 
dent of the direction of ki, rendering the two-point correlator of ^ in real space 
statistically isotropic. 

5. Non-FRW background expansion implies statistical anisotropy 

If the background metric is not FRW at some time, then the background metric 
is not FRW at t* because of assumption [TJ Thus (Gomez & Rodriguez, 2013), 

P^(ki) = P^^°{ki)[l + gs^{ki-ds^f+ gX{ki-d\f+ gl{ki-dl_)'^ + 

+ff+(ki • d+)^ + 5Ar(ki • Ni)^ -|- 5^(ki • N2)^ + 5Ar(ki ■ Ns)^], 

(25) 

where none of the g, except perhaps for the (if the background metric is again 
FRW at the time when is evaluated), are the same, and, in contrast, all the d 
directions are the same and equal to the preferred direction at t*. We conclude 
then that the power spectrum of ( does depend in this case on the direction of 
ki, rendering the two-point correlator of C in real space statistically anisotropic. 

6. Conclusions 

When we take sections and [ 5 ] together, they show that, under the estab¬ 
lished assumptions, the FRW background expansion is equivalent to the statisti¬ 
cal isotropy of C in connection with its two-point correlator in real space during 
some period of time that contains primordial inflation. 
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